黄土 高 原 粗 质地 土壤 剖面 水 分 运动 与 浅 层 地 下 水 补给 可 能 性 
模拟 
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摘 要 : 黄土 高 原水 土 流失 严重 ,生态 环境 脆弱 ,水 资源 短缺 ,地 下 水 对 保障 区 域 社会 经 济 发 展 和 维持 生态 系统 平衡 具有 
要 意义 ， 而 该 区 的 地 下 水 转化 和 补给 机 制 尚 不 明确 。 为 探究 黄土 高 原水 蚀 风 蚀 交错 区 土壤 剖面 深层 水 分 运动 及 降水 对 浅 层 
下 水 补给 的 可 能 性 , 利用 六 道 沟 小 流域 分 布 的 粗 质 地 风沙 土 样 地 2013 年 至 2016 年 土壤 剖面 0-600 cm 含水 量 数据 , 35 
HYDRUS-1D 模型 对 各 土 层 水 力 参 数 进行 反 演 和 验证 ， 并 用 于 模拟 样 地 土壤 深 剖 面 0-1 500 cm 水 分 运 移 过 程 。 结 果 显 示 , A 
平水 年 2014 年 (439 mm) 和 干旱 年 2015 年 (371 mm), 0-600 cm 土壤 含水 量 生长 季 末 与 生长 季 初 持平 或 略 有 亏 缺 ; 降水 充沛 久 
2013 年 (669 mm) 和 2016 年 (704 mm) 土 壤 含 水 量 生长 季 末 远 高 于 生长 季 初 ， 降 水 入 渗 深 度 超过 观测 深度 (600 cm)。 深 剖面 7 
分 运动 模拟 显示 , 2014 年 和 2015 年 剖面 含水 量变 化 不 明显 ,水 分 向 深层 运 移 微弱 缓慢 ; 但 是 , 2013 年 和 2016 年 降水 可 分 别 入 
渗 运 移 至 1 100 cm 和 1 200 cm 深度 ,， 远 超过 样 地 上 生长 的 旱 柳 根系 区 域 , 可 能 补给 浅 层 地 下 水 。 在 4 年 模拟 期 间 , 粗 质 地 风 
;> 十 水 分 入 渗 率 平均 为 606.28 cmd, 平均 土壤 蒸发 14.87 cm, 平均 植物 蒸腾 33.71 cm， 土 壤 中 的 水 分 主要 以 植物 蒸腾 形式 
损耗 。 在 2 个 丰 水 年 , 得 益 于 较 充 足 的 降水 和 粗 质 地 风沙 土壤 的 高 入 渗 率 ， 降 水 大 量 转 化 为 土壤 水 快速 向 下 入 渗 运 移 ， 模拟 
显示 当年 生长 季 末 降水 最 深 运 移 至 12m, 至 年 末 已 超过 模拟 深度 (15 m)， 水 分 继续 运 移 可 能 补给 浅 层 地 下 水 。 相 关 研 究 结果 
为 黄土 高 原水 蚀 风蚀 交错 区 地 下 水 来 源 和 补给 机 制 提 供 理论 依据 。 
关键 词 : 黄土 高 原 ; 水 蚀 风蚀 交错 区 ; HYDRUS-1D; 土壤 水 分 运 移 ; 数值 模拟 ; 深层 入 渗 
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A Simulation of water movement in coarse-textured soil and shallow groundwater recharge 
possibilities on the Loess Plateau, China“ 
Ge Jiamin!, Fan Jun!?**, Wang Sheng? 
(1.College of Resources and Environment, Northwest A & F University, Yangling 712100, China; 2.State Key Laboratory of Soil 
Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Chinese Academy of Sciences and 


Ministry of Water Resources, Yangling 712100, China) 


Abstract: The Loess Plateau suffers serious soil and water loss, resulting in a fragile ecological and environment and an intense 
deficit on available water resources. The groundwater is crucial to guarantee the regional social and economic development and 
maintain the ecosystem balance, however, the mechanism of groundwater transformation and recharge on the Loess Plateau is still 
not clear. The objective of this study was to investigate the water movement in deep soil profile and test the possibilities of rainfall in 
recharging the shallow groundwater in the water-wind erosion crisscross region on the Loess Plateau. Considering the restrictions of 
current equipment and techniques for deep soil profile moisture observation, HYDRUS-1D model was applied to calculate solute the 
soil hydraulic parameters and simulate deep profile 0-15 m soil water movement, by using of the data sets of meteorological factors 
and soil water contents in soil depth 0-600 cm in a Chinese willow field plot with aeolian sandy soil from 2013 to 2016. Upper 
boundary condition was set as atmosphere boundary conditions with inputs of rainfall and potential transpiration (Tp) and potential 
evaporation (Ep) which were partitioned from crops evapotranspiration (ETo), ETo derived from Penman-Monteith method and the 
partition into Tp and Ep had used the measured leaf area index (LAT) in different growing periods, and lower boundary condition was 
set free drainage. The optimized hydraulic parameters were derived from repeatedly inverse solution process, and then these 
parameters were used to simulate soil water movement in 0-1 500 cm, and the depth of soil water infiltration and the possibility of 
recharging shallow underground water were analyzed. Results showed that soil water contents at the end and the beginning of the 
growth season in 0-600 cm had a slight deficit in 2014 (439 mm, normal precipitation year) and 2015 (371 mm, dry year). However, 


there had a surplus of soil water at the end of the growth season as compared with the beginning of the growth season in 2013 and 


2016 (669 mm and 704 mm respectively, wet years), and the rainfall infiltration had exceeded the observing depth 600 cm. From the 
simulating results, soil water contents in deep soil profile had not change significantly in 2014 and 2015, and the water moves slowly 
into deeper soil layers, but rainfall could infiltrate to soil depths 11 m and 12 m, respectively, exceeding rooting zone at the end of the 
2013 and 2016, which indicated the possibilities of rainfall in recharging the shallow groundwater. During the simulated period of 4 
years, the mean value of water infiltration rate on aeolian sandy soil is 606.28 cm-d"!, average soil evaporation in 4 years is 14.87 cm, 
and average plant transpiration in 4 years is 33.71 cm, soil water contents were mainly lossed in the form of transpiration. Thanks to 
the high infiltration rate in coarse-textured soil in wet years, a large proportion of rainfall may transform into shallow underground 
water. Dry and wet years alternate on the Loess Plateau with intensive rainfall variations, the rainfall infiltration in wet years may be 
one way of shallow groundwater recharge. 


Keywords: Loess Plateau; Water-wind erosion crisscross; HYDRUS-1D; Soil water movement; Numerical modeling; Deep 
percolation 
地 下 水 补给 是 含水 层 和 含水 系统 从 外 界 获取 水 量 的 过 程 , 通常 受气 候 、 植 被 、 地 质 、 土 壤 、 地 貌 、 水 
位 埋 深 和 灌溉 等 下 垫 面 条 件 和 人 为 因素 的 综合 影响 , 因而 地 下 水 补给 过 程 十 分 复杂 。 黄 土 高 原水 蚀 风 蚀 交 
错 区 是 黄土 高 原 的 重要 组 成 部 分 ,面积 占 黄土 高 原 总 面积 的 28.56%, 水 蚀 风蚀 严重 , 水 资源 短缺 。 土 壤 质 
地 较 粗 ,是 黄河 下 游 河 床 粗 泥 沙 的 主要 产 区 趾 。 该 区 土壤 水 分 是 影响 区 域 生态 系统 生产 力 、 塑 造 区 域 植 被 
和 景观 的 重要 因素 外 。 在 当地 , 地 下 水 是 居民 生活 、 工 农业 生产 和 维持 区 域 生态 系统 平衡 的 重要 水 源 。 但 
T 目前 对 于 区 域 地 下 水 的 补给 和 转化 机 制 尚 不 明确 , 因此 研究 黄土 高 原 地 下 水 的 补给 对 水 资源 利用 、 地 下 水 
一 保护 和 模拟 水 文 过 程 有 重要 指导 作用 。 
CD 目前 对 深层 土壤 水 分 的 测定 和 保护 研究 很 多 , 不 同 土地 利用 模式 对 土壤 水 分 产生 较 大 影响 ,选择 合 
bs 的 植物 对 于 土壤 深层 水 的 保护 和 持续 利用 非常 重要 g。W. Hu 和 B. C. Si 在 时 间 变 异性 和 时 间 稳 定性 分 析 
- 基础 上 , 利用 SWC 模型 和 经 验 正 交 函数 ,估算 小 流域 尺度 上 空间 分 布 的 土壤 舍 水 量 。 李 蛤 等 外 利用 环境 示 
GN 踊 剂 气 估 测 长 武 城区 王 东 沟 流域 坡地 深层 渗 漏 量 , 平均 渗 漏 量 超过 20 mm a. 一些 学 者 在 巴 丹 吉林 沙漠 和 
M 日 古 乃 绿洲 利用 环境 同位 素 示 踪 手 段 , 根据 当地 气息、 土壤 质地 、 地 下 水 的 特征 及 水 化 学 条 件 , 研究 得 到 
= 地 下 水 补给 机 制 及 补给 来 源 , 并 定量 得 到 了 地 下 水 补给 深度 和 补给 量 53。 在 黄土 高 原 地 区 , 土 层 深厚 , 地 
Q 区 复 杂 。 诸 多 研究 显示 降水 是 黄土 区 地 下 水 补给 的 主要 来 源 &43。 但 降水 如 何 穿 过 深厚 的 黄土 层 补给 地 下 
水 的 过 程 尚 不 明确 中 I。 早 期 调查 指出 洛 川 扰 、 董 志 塘 等 黄土 城 区 主要 通过 裂 院 或 者 潜 刨 洞 补给 地 下 水 , 还 
在 裂 际 面 发 现 了 因 重 力 水 通过 形成 的 碳酸 盐 膜 而 较 难 通过 常规 入 渗 和 水 分 再 分 布 过 程 补给 04 刁 。 徐 学 选 、 
程 立 平等 0 指出 黄土 区 地 下 水 存在 大 孔 队 优先 流 补给 的 可 能 , 但 在 空间 上 并 不 普遍 ; 优先 流 的 发 生 与 土 
地 利用 方式 以 及 土壤 含水 量 有 关 , 土壤 干 层 会 减少 优先 流 的 发 生 ， 而 土壤 剖面 含水 量 较 高 时 则 容易 产生 优 
流 ， 从 而 对 地 下 水 进行 补给 &20。 区 域 “退耕 还 林 ( 草 ) 工 程 和 淤 地 坝 建 设 等 人 为 因素 导致 大 规模 植被 和 微 
地 形 改变 ,同时 引起 地 下 水 补给 和 转化 的 改变 09。 地 下 水 补给 量 在 不 同 区 域 不 同 研究 当中 存在 差异 ,例如 
黄土 高 原 无 定 河流 域 地 下 水 年 补给 量 为 11.38~15.69 mm， 占 多 年 平均 降水 量 (409 mm) 的 2.8%~3.9%07; 18 
JE ME 130 m 厚 层 黄 土 的 黑河 流域 地 下 水 年 补给 量 为 107+55 mm， 占 多 年 平均 降水 量 (571 mm) 的 
(19+10)%， 其 中 “优先 流 ” 占 (87+4)%08， 活 塞 流 由 于 干燥 层 的 存在 而 贡献 较 小 , 优先 流 为 主要 补充 机 制 。 
黄土 高 原 北 部 水 蚀 风 蚀 交 错 区 粗 质地 土壤 分 布 广泛 , 是 黄河 中 下 游 粗 泥 沙 的 主要 来 源 区 。 土 壤 孔 隙 度 
大 , 持 水 能 力 低 , 饱和 导 水 率 高 。 区 域 降水 量 在 108.6-891.1 mm 之 间 , 年 际 和 年 内 分 配 不 均 ,年 变 率 大 ， 降 
水 集中 且 暴 雨 多 ,两 强大 00。 有 研究 表明 ， 只 有 在 丰 水 年 降水 才能 垂直 入 渗 通 过 根系 活动 区 ， 从 而 对 地 下 水 
成 缓慢 补给 00， 大 的 降水 促进 了 大 孔 隐 优先 流 的 发 生 03。 入 渗 和 土壤 水 分 再 分 布 是 降水 、 地 表 水 、 土 壤 
水 和 地 下 水 相互 转化 的 一 个 重要 过 程 ， 研 究 年 际 不 同 降水 水 平 条 件 下 粗 质 地 土壤 深层 剖面 水 分 运动 , 对 揭 
示 降 水 对 浅 层 地 下 水 补给 可 能 性 、 水 蚀 风蚀 交错 区 地 下 水 补给 机 制 有 重要 意义 。 受 限于 水 分 观测 设备 和 实 
际 条 件 ， 更 深层 水 分 运动 很 难 实测 , 本 研究 采用 水 蚀 风 蚀 交 错 区 六 道 沟 小 流域 内 2013—2016 年 人 工 旱 柳 
(Salix matsudana) 林 地 土壤 0~600 cm 剂 面 含水 量 实测 数据 和 气象 资料 , 运用 HYDRUS-1D 反 演 获取 土壤 水 
力 参数 并 用 以 模拟 深层 土壤 水 分 运 移 ， 以 期 得 到 降水 在 粗 质地 土壤 中 入 渗 和 运 移 深 度 及 补给 浅 层 地 下 水 的 
可 能 性 , 为 揭示 水 蚀 风 蚀 交 错 区 土壤 剖面 水 分 运动 与 浅 层 地 下 水 补给 机 制 提供 参考 。 
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1 试验 材料 与 方法 
1.1 研究 区 概况 和 样 地 描述 

研究 区 位 于 陕西 省 榆林 市 神木 县 西 沟 乡 六 道 沟 小 流域 (110°21'~110°23 正 ，38°46'~38°51'N， 海拔 
1094-1274 m), 处 于 黄土 高 原 到 毛乌素 沙 地 的 过 渡 地 带 ， 以 及 流水 侵蚀 的 黄土 丘陵 区 向 干燥 剥蚀 作用 的 
哪 尔 多 斯 高 原 过 渡 的 水 蚀 风 蚀 交 错 带 ， 是 典型 的 生态 脆弱 区 。 该 区 属于 中 温带 半 干 旱 气 候 , 气候 变化 剧烈， 
年 平均 气温 8.4 C, 210 'C 活 动 积温 为 3 248.0 C, 无 霜 期 153 d, 全 年 日 照 时 数 2836 h, 年 总 太阳 辐射 量 
5 922 MJ.m-， 年 光合 有 效 辐射 2 900 MJ.m 习 ,年 平均 风速 为 2.5~2.7 ms-1。 多 年 平均 年 降水 量 437.4 mm, 年 
内 年 际 变 化 大 , 最 多 891.1 mm, 最 少 108.6 mm, 其 中 6 一 9 月份 的 降水 量 占 全 年 的 75% 以 上 , 年 平均 蒸发 力 
785.4 mm, 干燥 度 1.8, 属于 典型 的 半 干 旱 区 。 该 区 雨 热 同 季 , 植物 4 一 5 月 萌发 , 10 月 枯萎 。 
地 带 性 土壤 为 黑 坊 士 , 因 长 期 强烈 侵蚀 , 已 被 沙 黄 土 、 新 黄土 、 红 粘土 及 沙 地 上 发 育 起 来 的 风沙 土 、 
HEMER, 地貌 类 型 为 片 沙 覆 盖 的 梁 骨 状 黄土 丘陵 ,流域 内 地 形 非常 复杂 ， 沟 道 纵横 交错 ， 土 地 沙化 严 
重 。 退 耕 还 林 ( 草 ) 工 程 的 实施 , 使 该 流域 内 恢复 植被 大 量 增加 ,该 流域 目前 主要 人 工 或 自然 乔 灌木 有 杨 树 
(Populus simonii)、 蛙 柳 、 沙 柳 (Salix cheilophila) 和 柠 条 (Caragana Korshinskii), 主要 草 种 有 长 芒 草 (Stipa 
bungeana) 和 紫花 首 蒂 (Medicago sativa) 等 L191。 

试验 选择 旱 柳 样 地 为 研究 对 象 , 样 地 土壤 主要 为 粗 砂 质 土壤 , 在 当地 风沙 土 呈 片 状 、 斑 块 状 不 连续 分 
布 ,有 研究 表明 在 本 地 区 颗粒 组 成 大 于 0.05 mm E 69%P0， 也 有 研究 表明 0-200 cm 土 层 内 , 颗粒 组 成 小 于 
0.01 mm 占 13.8%PU。 在 风沙 土 丘 陵 区 , 0~200 cm 土 层 内 粗 砂 占 较 大 比例 。 据 实际 观测 , 0~260 cm 为 粗 砂 ， 
i 260-300 cm 为 粗 砂 细 砂 过 度 层 , 300-1 000 cm 之 间 820 cm 和 960 cm 附近 出 现 细 砂 与 黄土 混合 层 ， 其 余 为 
Lam 细 砂 质 土壤 , 据 此 认为 600-1 500 cm 土壤 为 细 砂 质 土壤 。 表 层 土壤 容重 为 1.62 g-cm?, 孔隙 度 为 38.996. 
~ 样 地 与 当地 气象 站 处 于 同一 海拔 高 度 1 216 m， 直 线 距 离 600 m。 样 地 上 的 旱 柳 是 当地 的 乡土 树种 ,树龄 近 
zh, 30 年 。 
~N 1.0 ”基础 数据 获取 

土壤 剖面 实测 含水 量 : 在 旱 柳 样 地 埋设 3 根 铝 制 中 子 仪 测 管 , 地 下 埋 深 均 为 600 cm, 采用 CNC503B 型 
中 子 仪 测量 土壤 剖面 体积 含水 量 (0, cms.cm-3), 在 2013—2016 年 4 月 至 10 月 每 隔 一 个 月 测定 一 次 。 因 土壤 
含水 量 浅 层 比 深层 变化 更 强烈 ,因此 在 0-100 cm 深度 土 层 , 测量 间隔 10 cm; 在 100-600 cm EE, 测量 间 
隔 为 20 cm. 

植被 根系 分 布 : 采用 80 mm 直径 根 钻 取 土 分 离 土 样 中 根系 ,获取 旱 柳 样 地 地 下 0~600 cm 深度 根系 分 布 
"一 特征 。 取 与 根系 吸水 密切 相关 的 小 于 2 mm 直径 细 根 生物 量 ， 分 析 土 壤 剖 面 根系 分 布 特征 ， 详 细 根系 取样 和 
F 处 理 过 程 见 前 期 研究 P2。 
^U UCBERUEHEZE BIET): 从 距离 样 地 直线 距离 600 m 的 试验 站 自动 气象 监测 站 获取 气象 数据 ， 主 要 
包括 1.5 m 高 度 净 辐射 量 (Rs,, MJ.m2h-D、 平 均 空 气温 度 (7，C)、 空 气相 对 湿度 (R 玉 , %)、2 m 高 度 风速 (Sw, 
m.s) 4 FE Jk Œ (P, mm) 55 . 3& Hj Ref-ET software (www.kimberly.uidaho.edu)， 选 择 内 置 ASCE 
Penman-Monteith RARS (ETo, mm.d-!)。 在 早 柳 不 同 生 长 时 期 用 SunScan 植物 冠 层 分 析 
仪 测 定 样 地 叶 面 积 指数 (LAD, 用 于 将 ETo 分 割 为 潜在 土壤 蒸发 (5E,,，mm-d-!) AUTE CE E ZEE, 
mm- dhA, 分 制 计算 公式 为 : 


d 


T, = ETo- SCF (1) 
Ep = ETo-(1- SCF) Q) 


SCF = 1 e" g "Pune LA (3) 


式 中 : SCF 为 植被 覆盖 系数 ,e 为 自然 指数 , rExtinct 为 冠 层 消光 系数 ,本 文 取 值 07505261, 


13 HYDRUS-1D 模型 
1.9.4 模型 介绍 


该 模型 包括 正解 模块 和 反 演 模块 ， 反 演 模块 允许 用 五 种 解析 表达 式 表 示 ] 
些 解 析 式 由 有 限 个 未 知 参数 组 成 。 模 型 模拟 可 以 方便 的 控制 环境 变量 和 边界 条 件 , 但 


ERKE Oh) K(h), 这 


般 需 要 复杂 的 参数 


输入 ,确定 合理 的 参数 比较 困难 。 反 演 过 程 中 先 给 这 些 参数 赋予 初 值 ， 然 后 把 赋予 初 值 的 gb 和 K(h) 代 入 


水 流 控 制 方程 ,并 结合 实验 的 初始 条 件 和 边界 条 件 ， 形 成 完整 描述 水 流 运动 状况 的 数学 模型 。 最 后 用 


Nash-Sutcliffe 效率 系数 Ens 和 均 方 根 误差 RMSE 评价 模型 的 模拟 精度 。 效率 系数 Ens 体现 了 模型 的 总 体 效 


K, 其 计算 公式 为 : 


D (P -O) 
n 一 2 
bo -0,) 


Ens=1- (4) 


式 中 : P; 为 预测 值 , 0; 为 观测 值 ，0, 为 观测 值 的 平均 值 ,n 为 观测 值 个 数 。Ens 的 变化 范围 -o~1。Ens 越 接近 


1， 则 模拟 效果 越 好 , 模拟 结果 可 信和 度 高 ; 通常 当 Ens>0.5 时 ， 即 认为 模拟 效 


果 可 接受 。 


RMSE 也 叫 回 归 系统 的 拟 合 标准 ， 表 示 的 是 模拟 值 与 实测 值 的 接近 程度 。 公 式 为 : 


RMSE = 


lY w(B-0y G 
nj 


RP: wi 表示 权重 。RMSE 越 小 ,代表 模拟 值 越 接近 观测 值 。 
1.8.0 ”初始 条 件 和 边界 条 件 


研究 区 沙 地 一 般 不 产生 径流 ,忽略 植被 截留 , 将 当年 降水 量 作为 水 分 输入 项 。 选 取 2013 年 至 2016 年 
实测 的 日 气象 数据 , 将 降水 数据 和 计算 得 到 的 E, 和 7 作为 大 气 边界 条 件 设 定 为 上 边界 条 件 ， 下 边界 条 件 设 


定 为 自由 排水 条 件 , 将 2013 年 0~600 cm 第 一 次 水 分 剖面 实测 值 作 为 初始 条 件 , 应 用 HYDRUS-1D 模型 进 
行 2 年 时 间 ( 共 730 qd) 的 水 分 运 移 模拟 ,首先 , 利用 反 解 模块 模拟 2013 一 2014 年 0-600 cm 剖面 水 分 变化 , 进 
行 参数 反 演 。 之 后 , 利用 反 解 获得 的 参数 , 将 2015 年 第 一 次 水 分 剖面 实测 值 作为 初始 条 件 , 模拟 得 到 的 
2015—2016 年 剖面 水 分 模拟 值 和 实测 值 进行 对 比 , 验证 参数 可 行 性 。 最 后 , 将 2013 年 0~600 cm 第 一 次 水 
分 剖面 实测 值 及 接近 600 cm 深度 的 含水 量 作为 600-1 500 cm 的 含水 量 作为 模型 初始 条 件 , 利用 验证 后 的 
参数 进行 2013—2016 年 0~1 500 cm 剖面 水 分 运动 模拟 .此 外 ,利用 观测 和 文献 数 志 


居 进 行 上 边界 通 量 验证 。 


根据 剖面 土壤 质地 差异 ,模拟 时 将 0-600 cm 土 层 剖面 划分 3 E, 分 别 为 0-260 cm. 260-300 cm. 


300-600 cm， 同 时 假设 600~1500 cm 土 层 土 壤 质 地 与 300~600 cm — $5 W} 
设置 模型 根系 分 布 参 数 , 不 考虑 根系 生长 。 
1.33 ”参数 确定 


5-3 IRA 


| 根 生物 量 分 布 规律 


采用 HYDRUS-1D 模型 的 反 演 模块 确定 剖面 土 层 水 力 参 数 : 进 气 吸力 的 倒数 (a)、 了 和 孔隙 大 小 分 布 系数 


(nm)、 饱 和 导 水 率 (Ks), 将 模拟 得 到 的 剖面 水 分 分 布 与 实测 值 进行 对 比 ,利用 


Nash-Sutcliffe 效率 系数 Ens 和 


均 方 根 误差 RMSE 检验 模拟 结果 和 实测 值 之 间 的 拟 合 程 度 ， 验证 获得 的 参数 和 入 渗 、 共 发 过 程 的 合理 性 , 不 
断 重复 反 演 过 程 ， 直 至 拟 合 效 果 最 优 且 参数 可 接受 , 得 到 最 优 的 反 演 参 数 。 之 后 将 确定 的 参数 代入 正解 模 


块 中 , 进行 2013—2016 年 的 深 剂 面 0~1500 cm 的 水 分 运动 模拟 。 

2 结果 与 分 析 

2.1 降水 特征 与 土壤 水 分 运动 
六 道 沟 流域 多 年 平均 年 降水 

试验 期 间 , 2013 年 和 2016 年 降水 


地 


ligi 


为 431.4 mm, 年 际 间 相 差 大 , 年 内 5 一 10 月 降水 量 约 占 全 年 降水 量 9296. 


分 别 为 669 mm 和 704 mm， 显著 高 于 多 年 平均 值 ， 为 丰 水 年 ; 2014 年 降 


水 量 439 mm, 接近 多 年 平均 值 , 为 典型 的 平水 年 ; 2015 年 降水 量 371.1 mm, 远 低 于 多 年 平均 值 ， 为 干旱 年 


(图 1)。4 年 内 最 大 日 降水 量 为 124.0 mm, 是 在 2016 年 7 月 8 日 ; 最 大 月 降水 量 也 在 2016 年 7 月 , 73 274.6 


mm。 由 于 试验 地 属于 粗 质地 的 沙土 , 土壤 水 分 含量 处 于 较 低 水 平 , 训 面 平均 含水 量 最 低 在 0.05 cm*cm?. 


土壤 水 分 年 内 受到 降水 的 强烈 影响 , 2014 E 2015 年 降水 量 少 ,土壤 剖面 水 分 含量 在 生长 季 初 和 生长 季 末 
整体 上 差异 不 大 , 说 明 当 季 降 水 均 已 通过 蒸 散发 或 者 渗 漏 等 损失 。 但 是 2013 年 和 2016 年 降水 较 多 , 土壤 
剖面 含水 量 在 生长 季 初 和 生长 季 末 有 显著 差异 ， 当 季 降 水 显著 提高 了 生长 季 末 剖面 大 部 分 层次 的 水 分 合 量 ， 
测量 的 最 深 处 生长 季 末 水 分 含量 仍然 高 于 生长 季 初 , 说 明 水 分 渗 漏 已 经 超过 了 试验 观测 深度 (图 2)。 
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1 2013—2016 年 六 道 沟 小 流域 月 降水 量 分 布 和 期 间 月 平均 降水 量 


Fig. 1 The monthly precipitation distribution and the monthly mean precipitation during the 2013-2016 


22 土壤 0-6 m 剖面 水 分 运动 模拟 

利用 2013—2014 年 反 演 得 到 的 参数 ， 模 拟 得 到 2015—2016 年 剖面 水 分 分 布 , 对 比分 析 2015—2016 年 
0~600 cm 剖面 水 分 生长 季 初 和 生长 季 末 实 测 值 和 模拟 值 。 结 果 显 示 剖 面 水 分 含量 模拟 值 与 实测 值 接近 , AE 
化 趋势 一 致 (图 2)。 同 时 分 析 了 50 cm、130 cm、300 cm、400 cm 和 540 cm EE 4 年 期 间 实 测 值 与 模拟 值 ， 
浅 层 土壤 受降 水 影响 较 大 , 在 50 cm 深度 土壤 水 分 随时 间 变 化 波动 较 大 ， 随 深度 的 增加 ,变化 幅度 逐渐 趋 
于 平缓 (图 3)。2013 一 2016 年 土壤 深度 模拟 结果 评价 系数 Ens 范围 在 0.439~0.739，80% 在 可 接受 范围 (Ens 
过 0.5) 内 , RMSE 在 0.002~0.009 cm? cm? WHA 1), 均 在 合理 范围 内 , 说 明 模 拟 结果 可 以 接受 。 

4 年 中 , 在 生长 季 初 和 生长 季 末 水 分 含量 在 0-280 cm 土 层 变化 明显 , 分 析 原 因 一 方面 可 能 是 受到 根系 
的 影响 ,根系 主要 分 布 在 0-100 cm 和 220~280 cm 土 层 , 分别 占 总 量 的 66% 和 2096071, d ZR PAS E dE 
含水 量 高 于 根系 多 的 土 层 。 另 一 方面 可 能 是 受 士 壤 剖面 质地 的 影响 , 280-300 cm 土 层 土壤 含水 量 显著 升 高 ， 
并 且 300 cm 土 层 以 下 土壤 含水 量 也 保持 在 一 个 较 高 的 水 平 (高 于 上 层 土壤 )。260~300 cm 是 粗 沙 向 细 沙 过 
度 层 , 很 有 可 能 使 得 水 分 下 渗 变 缓 而 造成 300 cm 左右 土 层 水 分 含量 增加 ， 并 且 柳 树 根系 主要 分 布 在 300 cm 
土 层 以 上 , 而 沙 地 大 量 的 水 分 入 渗 到 了 300 cm 土 层 以 下 , 这 说 明 柳 树 根 系 能 够 吸收 利用 的 水 分 较 少 ,超过 
根系 分 布 层 的 水 分 就 会 继续 下 渗 , 可 能 补给 浅 层 地 下 水 。 


土壤 体积 含水 量 土壤 体积 合 水 量 
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2 2013—2016 年 0-600 cm 剖面 水 分 实测 值 与 模拟 值 拟 合 程度 对 比 


Fig.2 Comparison between the observed and simulated soil volumetric water content profile in 2013—2016 
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图 3 2013—2016 年 不 同 土 层 实测 值 与 模拟 值 随 时 间 变 化 对 比 


Fig.3 Comparision between the observed and simulated soil volumetric water content at different soil depth in 2013—2016 


表 1 


率 定期 和 验证 期 不 同 土壤 深度 模拟 效果 评价 系数 Ens 和 RMSE 


Table 1 Goodness-of-fit assessment based on Ens and RMSE at different soil depth at the time of inverse and verification 


土壤 深度 率 定期 Inverse solution 验证 期 Verification period 
Soil layers (cm) Ens RMSE Ens RMSE 
50 0.652 0.006 0.731 0.009 
150 0.670 0.005 0.526 0.008 
300 0.439 0.006 0.613 0.007 
400 0.664 0.004 0.739 0.005 
540 0.483 0.002 0.576 0.004 


2.3 ”模型 验证 
HYDRUS-1D 模型 201 


范围 内 。 RMSE 在 0.006-0.016 cm?- cm? 范 


3 一 2016 年 不 同 测定 时 间 模 拟 评价 系数 Ens 值 从 0.433 到 0.915，63% 在 可 接受 


系数 说 明 模 型 模拟 0-600 cm 土壤 水 分 含量 结果 可 以 接受 。 


表 2 率 定 期 和 验证 期 不 同 测定 时 间 评价 系数 Ens 和 RMSE 


Table 2 Ens and RMSE evaluation coefficients at different measured time of correction and verification 


围 内 , 均值 为 0.009 cm?- cm? CX 2), 2013 年 至 2016 年 的 两 个 评价 


期 (年 -月 -日 ) 率 定 期 Correction 日 期 (年 -月 -日 ) 验证 期 Verification 
Date (year-month-day) Ens RMSE Date (year-month-day) Ens RMSE 
2013 年 4 月 24 0.433 0.007 2015 年 4 月 21 H 0.644 0.007 
2013 4£ 6 H 14 E 0.586 0.008 2015 Æ 6 H 23 0.82 0.006 
2013 4E 8 H 16 0.625 0.007 2015 £7 H 20 H 0.468 0.012 
2013 4 10 H 15 H 0.451 0.016 2015 年 10 H 23 0.455 0.008 
201445 H 15 0.574 0.012 2016 ££ 4 H 22 H 0.516 0.008 
20144 6 H 15 H 0.652 0.011 2016 £5 H 16 0.463 0.009 
2014 4 8 H 20 0.552 0.010 2016 4E 8 H 25 0.696 0.010 
2014 4£ 9 H 19 0.477 0.012 2016 £ 9 H 23 H 0.915 0.007 


模型 输出 的 2013—2016 年 结果 ( 表 3) 显 示 ， 


2016 年 的 丰 水 年 , 生长 季 末 
大 于 0, 与 表 2 对 应 。 


同时 , 经 模型 反 演 得 到 的 土壤 水 力 参 数 可 以 月 


土壤 剖面 含水 量 明 


T 


表 3 2013—2016 年 模型 输出 地 表 水 通 量 


Table 3 The surface water flux outputs of the simulation model in 2013—2016 


条 水 量 与 实际 蒸 散 量 之 差 与 实际 情况 一 致 ， 即 在 2013 年 和 


显 比 生 长 季 初 剖面 含水 量 增 加 ， 降 水 量 与 实际 花 散 量 之 差 应 


来 模拟 更 深层 次 土壤 水 分 运动 ( 表 4). 


降水 量 - 蒸 散 量 


年 份 土壤 蒸发 植物 蒸腾 PEHE ZSE - ME 植物 蒸腾 / 蒸 散 
: mE M 降水 量 Precipitation minus 
Evaporations Transpirations — Evapotranspirations MM A" Transporates 
Precipitation(cm)  evapotranspirations D 
Year (cm) (cm) (cm) /Evapotranspirations 
(cm) 
2013 13.065 46.132 59.197 66.94 7.143 0.779 
2014 17.548 30.338 47.886 43.92 —3.966 0.634 
2015 15.038 22.599 37.637 37.11 —0.527 0.600 
2016 13.810 35.711 49.581 70.43 18.429 0.721 


表 4 模型 在 率 定期 间 输 出 的 土壤 水 力 特征 参数 


Table 4 Soil hydraulic characteristic parameters from the model simulation during the period of inverse correction 


剖面 层次 凋 蓉 含 水 0r 饱和 含水 量 进 气 吸力 倒数 — 孔隙 分 布 系 饱和 导 水 率 Ks 
Soil layers(cm) (cm?- cm?) 0 s(cm?-cm?) a (cm!) 数 (n) (cm:d-!) 
0~260 0.025 0.399 0.043 3.117 712.8 
260~300 0.019 0.371 0.038 1.783 629.8 
300~600 0.015 0.366 0.032 1.697 476.2 


2.4 ”土壤 深 剖 面 0~15 m 水 分 运动 模拟 

模型 通过 验证 之 后 ,保持 相关 参数 不 变 , 将 模拟 深度 增加 至 1 500 cm， 模 拟 土壤 水 分 向 深层 的 运动 过 
程 (图 4)。 

2013 年 降水 量 较 高 , 模拟 结果 显示 ， 雨季 末 的 10 月 份 300~800 cm 的 剖面 水 分 
加 , 10 月 份 入 渗 到 1 100 cm 左右 ， 因 为 生长 季 初 (4 月) 与 生长 季 末 (10 月 ) 剖 面 水 分 含量 分 布线 在 1 100 cm 重 
E, 水 分 含量 趋 于 0.08 cm*cm?. HF 10 月 份 之 后 , 继续 进行 模拟 , 结果 显示 , 300 cm 以 下 土壤 水 分 继续 
缓慢 下 渗 , 100 d 之 后 ,水 分 入 渗 深 度 超过 
1 500 cm. 

2014 年 季 初 、 季 末 及 年 底 3 个 时 间 点 剖面 水 分 分 布 变化 不 大 , 但 是 随时 间 延 长 , 深层 剖面 水 分 含水 量 
分 布线 没有 变化 , 可 能 是 2013 年 较为 丰沛 的 降水 在 深层 继续 下 渗 ， 而 浅 层 含水 量 分 布线 在 400 cm 左右 就 
重合 , 说 明 当 季 降 水 量 并 没有 下 渗 超 过 400 cm. 
2015 年 属于 干旱 年 , 在 前 一 年 平水 年 的 基础 上 ,土壤 剖 面 水 分 含量 处 于 较 低 水 平 。 从 模拟 结果 看 , 没 
有 明显 的 土壤 水 分 下 渗 , 因此 2015 年 降水 无 法 补给 到 浅 层 地 下 水 。 但 是 从 2014 FERE 2015 年 季 初 的 4 
个 多 月 期 间 , 300 cm 以 下 土壤 剖面 水 分 含量 整体 有 降低 趋势 ,这 说 明 2013 年 的 降水 , 经 过 2014 全 年 和 2015 
年 的 生长 季 前 这 段 时 间 的 入 渗 再 分 布 ， 土壤 水 分 的 入 渗 深 度 已 经 超过 1 500 cm, 水 分 可 以 继续 入 渗 , 最 终 
补给 浅 层 地 下 水 。 

2016 年 是 降水 量 较 高 的 丰 水 年 , 10 月 份 300~700 em 剖面 水 分 含量 相 比 4 月 份 明 显 增加 , FRAK 

$ 量 分 布线 在 1200 cm 左右 与 季 初 水 分 含量 分 布线 重合 ， 即 季 末 水 分 入 渗 到 
1 200 em, 模拟 到 年 末 水 分 含量 分 布线 延伸 超过 1.500 ecm。 按照 2013—2015 年 期 间 的 水 分 运动 状况 ， 本 年 


度 的 土壤 水 分 将 在 2017 年 和 2018 年 生长 季 之 前 ， 有 更 多 的 水 分 入 渗 超 过 1 500 cm 深度 ,向 更 深层 运动 。 
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Fig.4 Soil water content at depth 1 500 cm at the beginning of the growing season, the end of the growing season and the end of the year in 
2013-2016 
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3 讨论 

根据 前 期 研究 , 该 样 地 2013 年 5 一 9 月 土壤 蒸发 总 量 (5 2)73 17.0 em, 植物 蒸腾 总 量 (7 为 26.2 cm", T 
iJ ET AS 为 0.606, 2015 年 Tg/ET gs 为 0.369281。 模 型 模拟 结果 得 到 2013 年 、2014 年 、2015 年 和 2016 EHI TIET 
分 别 为 0.779、0.634、0.600 和 0.721, 4 年 的 TVET 平 均值 为 0.684, 模型 模拟 值 与 已 有 测定 值 接近 , 说明 模 
型 模拟 的 上 边界 通 量 结果 与 实测 接近 。 降 水 进入 土壤 后 , 除了 上 边界 的 蒸 散 发 损失 之 外 ,就 是 在 土壤 中 的 
储存 和 深层 渗 漏 。 雨 强大 ,历时 短 , 一 次 性 降水 多 是 黄土 高 原 典 型 的 降水 特点 C9。 对 于 有 高 孔隙 度 的 砂 质 
土壤 来 说 ,， 其 饱和 导 水 率 高 ,田间 持 水 量 低 ， 强 降水 后 水 分 入 渗 迅 速 。 降 水 结束 后 ， 强 烈 燕 发 使 砂 质 土壤 表 
面 很 快 形成 干 沙 层 B0,， 大 孔 隐 毛管 上 升水 通道 断裂 ， 供 水 作用 停止 ， 随 后 土壤 蒸发 仅 发 生 在 深层 土壤 , 水 
汽 通过 土壤 孔隙 再 扩散 到 大 气 中 去 , 茹 发 速率 比 同样 条 件 下 水 面 蒸发 小 得 多 BY, 干 沙 层 减少 了 地 下 水 受 强 
烈 蒸发 的 损失 ,起 到 了 保护 深层 土壤 水 分 的 作用 。 当 入 渗 深 度 超 过 根系 深度 时 ,土壤 水 分 不 受 蔡 散 发 影响 ， 
继续 向 下 渗 漏 0021。 在 这 种 条 件 下 ,即使 是 干旱 年 , 也 有 可 能 补给 浅 层 地 下 水 。 

在 丰 水 年 ， 由 于 土壤 质地 原因 , 大 的 降水 很 难 在 根 区 滞留 而 被 植物 利用 ， 尽 管 一 定时 段 内 降水 量 可 能 
小 于 同期 植物 需 水 量 , 但 仍 有 水 分 通过 根 区 下 渗 至 深层 土壤 00， 否 则 降水 转化 成 士 壤 水 后 不 参与 地 下 水 循 
H, 通过 蒸腾 燕 发 作用 返还 到 大 气 中 B2。 模 拟 结果 显示 ,2013 年 和 2016 年 降水 量 相对 薰 散 量 有 乔 余 ， 目 
2016 年 降水 量 显著 高 于 同期 蒸 散 量 ， 入 渗 深 度 也 超出 了 根 区 范围 , 通过 降水 入 渗 与 水 分 再 分 布 过 程 ， 降 水 
可 能 补给 到 浅 层 地 下 水 , 进入 小 流域 地 下 水 循环 的 过 程 中 。 

Mathieu 和 BariacB3] 认 为 ， 通 过 对 比 降 水 和 土壤 水 中 同位 素 变化 特征 , 将 降水 入 渗 分 为 2 个 过 程 : 一 是 
降水 入 渗 到 土壤 基质 中 , 与 土壤 浅 层 自 由 水 完全 混合 后 继续 下 渗 ,， HARR E RAAM, 二 是 降水 
bom, 通过 土壤 大 孔 际 以 优先 流 的 形式 快速 通过 浅 层 土壤 ， 最终 导 致 深层 渗 漏 .Gazis 和 FengB4 的 研究 显示 , (ETE 
^ 层 土 壤 ， 降 水 在 入 渗 过 程 中 以 活塞 流 形式 下 渗 ， 同时 部 分 雨水 能 够 以 优先 流 的 形式 快速 抵达 深层 土壤 ， 这 
aia 种 优先 流 的 存在 , 使 得 降水 能 够 快速 补给 浅 层 地 下 水 。 但 最 近 关 于 地 下 水 的 同位 素 研究 表明 ， 当 季 降 水 可 

能 并 不 容易 和 浅 层 的 土壤 水 分 充分 混合 ,， 并 以 活塞 流 的 方式 下 渗 , 而 是 以 “ 易 动 水 ”的 形式 直接 补给 地 下 水 ; 
只 有 原来 保持 在 土壤 根 际 附近 的 小 孔隙 中 的 “ 难 动 水 ”被 植物 耗竭, 根 区 才 有 可 能 得 到 当 季 降水 的 补给 B4。 
目前 , 大 多 数 研究 认为 黄土 高 原 的 震 区 丰富 的 地 下 水 是 通过 优先 流 补给 的 , 因为 传统 “活塞 ” 流 穿 过 厚 层 
黄土 需要 几 十 甚至 几 百 年 的 时 间 凤 3539。 而 在 丘陵 沟 殖 区 ， 降 水 可 以 通过 复杂 地 形 形 成 的 特殊 通道 补给 浅 
层 地 下 水 ,例如 黄土 和 十 土 壤 的 接触 面 , 边 坡 裂隙 等 81。Yasuda 等 63 在 六 道 沟 流域 的 研究 表明 ， 降 水 以 坡 
e 面 径流 的 形式 在 淤 地 坝 沟 头 部 位 汇集 入 渗 , 进而 补给 浅 层 地 下 水 。 
= 与 试验 样 地 临近 的 沟 道 基 流 观测 结果 显示 , 2015 年 6 一 10 月 的 基 流 流量 随 降 水 的 变化 而 变化 (图 4), 但 
《9 是 基 流 的 增加 对 降水 响应 存在 滞后 现象 ,这 与 Tan 等 67 研 究 结论 一 致 。 如 果 6 月 28 日 .7 月 17 日 .8 月 3 

日 、9 月 4 日 、9 月 8 日 、10 月 2 日 的 基 流 曲线 峰值 分 别 是 对 8 月 28 日 9 月 10 日 9 月 28 日 、10 月 19 
日 、10 月 25 日 、11 月 5 日 的 降水 的 响应 , 那么 间隔 时 间 分 别 为 61 d. 55d. 59d. 45d. 47d. 33d, 平均 
滞后 时 间 大 约 为 50 d。 基 流 可 能 来 源 处 的 地 表 与 基 流 出 流 处 岩层 相对 高 差 为 17 m, 计算 得 到 土壤 水 分 运动 
平均 速率 为 34 cm d-!。 试 验 观 测 的 旱 柳 样 地 距离 基 岩 30 m 左右 , 按照 这 个 速率 ， 季 末 降 水 期 基本 结束 后 ， 
将 近 3 个 月 的 时 间 ， 降水 可 以 到 达 基 岩 , 进而 可 能 补给 到 浅 层 地 下 水 。 土 壤 剖 面 水 分 运动 模拟 结果 显示 ， 雨 
季 结 束 后 , 2013 年 和 2016 年 两 个 丰 水 年 份 , 水 分 下 渗 至 基 岩 处 需要 100 d 左右 的 时 间 , 二 者 的 入 渗 速 率 基 
本 接近 , 证 明 降 水 通过 土壤 基质 入 渗 与 再 分 布 过 程 也 可 以 补给 该 区 的 浅 层 地 下 水 , 因为 HUYDRUS-1D 模 
拟 的 水 分 运动 过 程 并 未 考虑 优先 流 现象 。 但 是 , 在 水 蚀 风 蚀 交 错 区 ,黄土 与 风沙 土 交 错 分 布 的 地 貌 特征 下 ， 
风沙 土 导 水 率 高 ， 持 水 性 差 , 区域 上 成 为 降水 补给 地 下 水 的 “优先 通道 ”。 在 风蚀 水 蚀 交 错 区 , 往往 在 风沙 
士 堆 积 的 附近 沟 道 有 “泉水 ”， 就 是 例证 。 
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5 2015 年 基 流 与 日 降水 变化 
Fig.5 Variety of base flow and daily precipitation in 2015 
4 结论 
黄土 高 原水 蚀 风蚀 交错 区 降水 年 际 差异 大 , 在 连续 4 年 观测 期 间 ， 最 大 年 降水 量 约 是 最 小 年 降水 量 的 2 
党 ， 水 分 入 渗 过 程 受降 水 量 影 响 明 显 。 在 区 域 平水 年 或 干 早年 , 粗 质地 土壤 剖面 含水 量 仅 在 3 m 以 内 深度 波 
动 , 3 m~6 m 水 分 含量 变化 微弱 。 但 是 , 丰 水 年 2013 年 和 2016 年 生长 季 末 训 面 水 分 含量 显著 高 于 生长 季 初 


得 益 于 较 充 足 的 降水 和 粗 质地 土壤 的 高 入 渗 率 , 降水 大 量 转化 为 土壤 水 快速 向 下 入 渗 运 移 .在 观测 期 间 的 2 
个 丰 水 年 ， 降 水 入 渗 深 度 超过 观测 深度 (6 m)， 模 型 模拟 显示 当年 植物 生长 季 未 降水 最 深入 渗 运 移 至 12 m, 
年 未 已 经 超过 模拟 深度 15 m， 远 超过 根 层 , 水 分 继续 向 深层 运 移 ， 可 能 补给 浅 层 地 下 水 。 该 区 黄 十 与 风沙 十 
交互 堆积 形成 的 特殊 地 貌 条 件 下 ,风沙 土 剖 面 可 能 是 浅 层 地 下 水 补给 的 “特殊 通道 ”。 
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